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A scheme for calculating corrections to all orders to the entropy of any thermodynamic system 
due to statistical fluctuations around equilibrium has been developed. It is then applied to the BTZ 
black hole, AdS Schwarzschild black Hole and Schwarzschild black Hole in a cavity. The scheme 
that we present is a model independent scheme and hence universally applicable to all classical black 
holes with positive specific heat. It has been seen earlier that the microcanonical entropy of a system 
can be more accurately reproduced by considering a logarithmic correction to the canonical entropy 
function. The higher order corrections will be a step further in calculating the microcanonical 
■ entropy of a black hole. 

o : 
<n : 

\ I. INTRODUCTION 

<: 

. It is well known that black holes obey laws which are analogous to the laws of thermodynamics 0. The similarity 

£C) 1 of the expression 

dM = T H ^- + $dQ (1) 

, with the expression in thermodynamics 

o ■ 

O 1 dE = TdS - PdV (2) 

and the property of the nondecreasing black Hole area found by Hawking^ strongly supported the Bekenstein 
conjecture [lj that the entropy of a black hole is proportional to its area. The event horizon prevents us from knowing 
O any other information except the mass, the charge and the angular momentum of the black hole (commonly called as 
the "No hair theorem" 3]). Thus various black holes formed by different configurations of infalling matter can have 
!— i . the same mass, charge or angular momentum. These variables are like the thermodynamic variables of energy and 
pressure. In thermodynamics we have to trade information to describe the system's overall behaviour. A system can 
, have various possible configurations leading to the same overall behaviour. This leads to the concept of entropy. It is 
T ■ proportional to the logarithm of the number of possible configurations. 

Various attempts have been made to explain the microscopic origin of the entropy of the black hole in terms of 
such degrees of freedom All of them use the Boltzmann definition of entropy 

S = lnfl{E) (3) 

where £l(E) denotes the number of different configurations a system can have and Ub, the Boltzmann's constant has 
been set to 1. The degrees of freedom for a black hole are hidden from us due to the event horizon and their logarithm 
can be interpreted as the entropy of the black hole, in the same spirit as the entropy of a thermodynamic system. 

In statistical mechanics we have the notions of the micro-canonical and the canonical ensemble. In general one 
can define entropy in both these ensembles. An important thing that should be noted is the difference between the 
two definitions. While in the microcanonical ensemble the energy is fixed at a given point E, in the canonical ensemble 
the energy is allowed to fluctuate about a mean energy (E) . If one uses the notion of entropy as information, one can 
conclude that the canonical entropy should be more than the microcanonical entropy. This is because in the canonical 
ensemble there is an additional ambiguity which remains about which states a system can take as the system can be 
in a state with energy close to (E) apart from the states with energy equal to (E) . 
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In the next section we show that the leading order difference between the canonical and the microcanonical 
entropy for any thermodynamic system is given by — ^ ln(CT 2 ) where C is the specific heat of the system and T is 
its temperature as has been shown earlier This is accounted by the logarithmic corrections found out by other 
methods |14| suggesting that the semiclassical arguments leading to Bekenstein Hawking entropy give us the canonical 
entropy That such a correction to the canonical entropy leads closer to the microcanonical entropy also has been 
verified analytically and numerically in [Tsl ]. 

It was conjectured in (^^| and references therein) that the next order correction should be proportional to the 
inverse of area 19]. We develop a scheme to calculate corrections to all orders to the Bekenstein Hawking entropy 
using techniques of thermodynamics. 

The rest of the paper is organised as follows. In section ITTT1 we apply our results to the BTZ blackhole, the Anti 
de-Sitter Schwarzschild black hole and the Schwarzschild black hole in a cavity. In section lTvl wc summarise the results 
obtained by us, explain the various assumptions which went in and the future directions of work. 



II. ANALYSIS OF THE THERMODYNAMIC FLUCTUATIONS 



Consider a many body system with single particle energy spectrum E n completely specified by the quantum 
number n 20]. Thus we have 

E n = f(n), n e Z* (4) 

where we assume that / is a bijection from IN to R . We also assume that f(n) has a differentiable inverse F such that 
n = F(E). Let the degeneracy of the states at energy E n be given by £l(E n ). The quantum density of the system is 
then given by 

p(E) = J2n(E n )5(E-E n ) . (5) 

n 

We use the delta function identity 

5{E - E n ) = 5{E - /(n)) = 8{n - F{E))\F'(E)\ , (6) 

to transfer the delta function to n from E. Here the prime denotes the first derivative of F with respect to E. The 
quantum density of states is then given by 

n— oo 

p(E)=n(E)\F'(E)\ "£6(n-F(E)). (7) 

n=0 

The Poisson summation formula 

oo oo 

J2 5{n - F(E)) = exp(2iirkF(E)) . (8) 

n— — oo k— — oo 

leads to 

oo oo 

^T5(n- F{Ej) = 1 + 2^ cos{2irkF{Ej) . (9) 

n=0 fc=l 

considering the fact that F(E) > 0. Using Eq.0 and the Poisson summation formula we get 



p(E) = Q(E)\F'(E)\ ( 1 + 2^cos(2tt£;.F(£0) ) 



(10) 



\ k=i 

Note that the notation and the derivation till here follows the scheme of [2l|. We thus obtain 

p(E) = Sl(E)\F>(E)\ , (11) 

where p{E) denotes the averaged smoothed density of states. Now we consider the case of N noninteracting particles 
constituting the many body system. The microcanonical entropy is just then the natural logarithm of the degeneracy 
D,{E n ) which is 

S = Hp(E)\F'(E)\- 1 } , (12) 
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If one assumes that the function p(E) is smooth and neglects the oscillating components which arise from the 
discreteness of the energy levels (which can be done for a macroscopic system with large energy E), then 

p(E) = p(E) (13) 

We now consider the canonical ensemble to calculate the density of states. The canonical N particle partition function 
is given by 

/>oo 

Zc{(3) = exp(-/3£n) = / p{E) exp(~f3E)dE , (14) 

where p(E) is defined by Eq.JSJ. The canonical ensemble is used as a tool here to extract the canonical density of 
states 22] . We observe that the canonical partition function is the Laplace transform of the density of states. So we 
obtain p{E) by Laplace inversion with respect to /3|23t,l24j. which will be our variable along the imaginary axis. Thus 



p{E) = -L / Z c (f3)e^ E dp (15) 



2iri 

which can be rewritten as 

-i pioo 

p{E) = — / J»{Zc{fi))+f>E d p (16) 
2tt« J_ ioo 

Here E is the canonical ensemble energy average E. We expect that the function S{fS) = ln(Z c (P)) + f3E has 
an extremum which the thermodynamic system will tend to. Let the extremum be at (3 = /3q. One can expand the 
function S{0) in a Taylor series around /3q as 

S(0)=ln(ZM) + 0E = S o (E) + ± (J^\nZ c ^ [fi - /3 ) 2 + ± (J^ lnZ c ) (/3 - /3 ) 3 + .... (17) 

Since the first derivative of S with respect to vanishes at /3q using Ea. (|16fl and Ea. (|17f) . we get 

pioo 

Note that here we have assumed that the second derivative is positive and defined as : 

d 2 



p (E) = — I e So(E) e ± a t(f3-M 2 + ±a 3 (f3-Pof + { ! a i (f3-f3 ) i + ... ^g) 



as) 



while for the rest n, a n denotes the n-th derivative of In Z c 

_ ( d n 

an - \W 

Changing variable from (/? — /3q) to iy, we get 



lnZ c \ (20) 



pi E) = ±-e s °W [ e-^^ 2 e^3^ (21) 



— oo 



The second exponential can be expanded as a series noting that for a large thermodynamic system the exponent is 
small to give 



2tt I \ ' m! , 

\m=0 \n=3 
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As we know that the L.H.S is a real quantity we can safely ignore the terms giving imaginary quantities in the 
R.H.S. The complete expansion will involve the multinomial expansion. We use terms only up to second order in m 
here i.e up to m = 2, noting that the highest order in the neglected next terms is 0(y 12 ) on account of considering 
only non-imaginary terms. (The cube of the term with n — 3 does not contribute being imaginary.) To simplify matters 
just consider the expansion which we did in the last step. 



ml ( E/ 



\ n — 3 




(23) 



The terms which contribute real terms are then (up to m = 2) 



E 



a 2 n(-l)V 
(2n)! 



oo oo 



EE 



a n a m (i) n+m y n+m 

71 ' 771 I 



n— 3 rn — 3 



n -f m = 2k k e IN 



(24) 



Using Eq.J22J) and the standard integrals 



e~i a 2V 2 dy = 



(25) 



and 



V zn dy = 



/27T (2n- 1)!! 



v 2n 



(26) 



we get the integrated expression for p(E) as follows 



p(£) 



2tt 



E 



a2n(-l)" 

(2n)!!af" 



OO OO 

EE 



a„a m (-l) fe (2A:-l)!! 



2! 



n— 3 m — 3 



(27) 



The above derivation is in the spirit of corrections to the saddle point formula 25, 26]. Here we would like to contrast 
our approach with the approach in [27| where the /3 3 term has been ignored as it contributes an imaginary term. 
Hence some extra contributions are missing following the approach in [2jJ. 

In the limit of a macroscopic system with a very large value of energy, we neglect the oscillating components in 
p(E) and put it equal to p{E). To get the microcanonical canonical entropy we take the logarithm of both sides to 
get 



Hp{E)) = In(n(£7)) + HF{E)) = S{E) + ln(F(£?)) 



(28) 



S(E) = So(E) - HF(E)) - i Hal) + In 1 + £ + 4 E E +oft (29) 

Next we use the approximation 

ln(l + x) ~ x (30) 

to get the final expression 



even cm , 1 w „2n , y> Q 2 »(-l)" 1 ^ ~ a„q m (-l) fc (2fc - 1)!! 

= So(£) - In(F(£0) - 5 M°a) + E ^n^T + 5? E E + O(^) (31) 



n— 3 m— 3 



This is our formula to calculate the corrections to black hole entropy which is valid for any thermodynamic system. 
Thus in addition to the logarithmic correction the microcanonical entropy formula has, a logarithmic correction 
dependent on the energy spectrum of the black hole. Such corrections and their possible origins have been studied in 
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III. APPLICATIONS TO BLACK HOLE ENTROPY 

We now apply Eq.JJlJ to find out entropy corrections to the BTZ blackhole, the Anti de-Sitter schwarzschild black 
hole and the Schwarzschild black hole in a cavity. 

A. BTZ Black Hole 

Consider a non-rotating BTZ black hole in three dimensions with metric 

ds 2 8M^j dt 2 +(j2~ 8M^j dr 2 + r 2 d9 2 . (32) 

The black hole behaves as a thermodynamic system with the thermodynamic quantities namely its Hawking 
temperature and Bekenstein-Hawking entropy given by the expressions 



T = I±- 

H 2tt£ 2 



1 



ir 2 £ 2 



So (33) 



So= 2 -^ (34) 





- 2 JL 


9 ^0 

" K 2 


Kn\ 


- (-i)" 


n!S*™ +1 


Ph +1 ' 


K n 



where the horizon radius is denoted by r + and given by r + — y/SG^Mi, M being the mass of the black hole and 
I is related to the cosmological constant by A = — i. It can be shown that M, Th and S satisfy the first law of 
Thermodynamics. 
Hence we have 

S ((3)=^f ± =* (35) 

PH PH 

where K is a constant. 

Using this expression we can calculate the required a„'s. 

TS C3 

(36) 

K)n>a = (-1)"^t - i-lT'^- (37) 
Ph 

Thus using Eq. i|31|) we find that in addition to the logarithmic correction arising due to the energy spectrum we have 
another logarithmic correction given by 

-\Ma 2 ) = -hnS n + K 1 (38) 

where Ki is a constant. This reproduces the correction in |l8j . 
Now we proceed towards the higher order terms : 

Higher order corrections = £ + ^yEE ~ ~ nXmXaf ( 39 ) 

n=2 v ' 1 n=3m=3 2 

We consider these terms one by one and use Ea . I|36|) and Ea. 137(1 to get 

f a* n (-ir = f My (2n- 1)11 

2_, (2n ) !!a 2n 2^\2j S"- 1 [ j 

a n a m (-l) k (2k 1)!! ^ ^ (-l) fc (2k 1)!! 
2_, n \m\af ^ ^ 2 fc +i S k ~ 2 1 ' 

Thus we have corrections as polynomials in -J-0|"l23- The leading term in them being appearing with the 
coefficient occurring for n=2 for the first of the equations above and k — 3 translating to m — 3, n — 3 for the 
following equation. Thus we get 

S(E) = S (E) - \ ln(5„) - + ... - HF(E)- 1 ) (42) 
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Also note that we can trust the factor of the S 1 term as the 0(y 12 ) terms neglected in the expansion of the 
exponential in Eg. 112411 give a term -i 



B. Anti de-Sitter Schwarzschild Black Hole 

Consider a d-dimcnsional Anti de-Sitter Schwarzschild black hole. We parameterise the cosmological constant in 
d-dimensions in terms of the length scale following as 

(d-l)(d-2) , 
A = - W 2 " ^ 



The AdS Schwarzschild black hole has the metric 

WttM r 2 \ , , ( 16tt 



<l»- = - { 1 : — r— + — i dt 2 + ( 1 - + ) dr 2 + r 2 dn 2 d , ( 44 ! 

V (d - 2)n d - 2 r d - 3 I 2 ) \ (d-2)Oi_ 2 r d - 3 I 2 ) d ~ 2 V ; 

The parameter M is the conserved charge associated with time-like Killing vector. In terms of the horizon radius 
then the entropy and the Hawking temperature are given by 

So = 4 + (45) 



T H = 



(d - l)r% + (d - 3)1 
4:Trl 2 r + 



(d - 3) n d _ 2 



4tt 



1 



(d-1) 



(d - i)i 2 



S (d " 2) ee fc 1 - i 1 S (d - 2) (46) 



where fci is a constant. The last approximation arises considering the large r+ limit so that r + 3> I which is the 
so-called high temperature limit. The specific heat then turns out to be |l8j | positive for a wisely chosen value of I. 
In the limit of I — > we get the specific heat to be positive for d > 3 and given by the expression 

C = (d - 2)S (47) 

The relation between So and (3h is then given by 

So = k d ~ 2 [3 H {d ~ 2) = k 2 (3 H {d - 2) (48) 

where k 2 is a constant. Using this we calculate the corrections to the entropy. For the first order corrections we need 
the second derivative of the entropy with respect to inverse temperature. In our notation 

a 2 2 = k 2 {d-2)(d-l)P H d (49) 

a 3 = k 2 (d-2)(d-l)(-d)f3 H (d+1) (50) 

a 4 = k 2 (d - 2)(d - l)(d)(d + l)f3 H (d+2) (51) 
And with the usual notations that has been followed up till now we have, 

% oc So 1 - (52) 

^Sr « s 2 - k (53) 

the neglected term 0(y 12 ) in the exponential expansion goes here as 0(4j) again as can be verified and thus is less 
than the one we are talking about. Thus the corrected entropy up to order ^- is 

S = S '-«-VA,„ 4+ _fi_^_^ (54) 

For example for a four dimensional Anti de-Sitter Black hole (d=4) we find that 

S = So -MF(E)- 1 )- hi S + (55) 

3d Sq 
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C. Schwarzschild in a cavity 



A Schwarzschild black hole has a negative specific heat but a thermodynamically stable Schwarzschild black hole 
can be obtained within a finite cavity by immersing it in a thermodynamic bath|34|. Let us find what are the higher 
order corrections to the entropy of such a black hole. 

Again we proceed the same way as the cases above. For such a black hole we have 

S = AttM 2 (56) 



where R is the radius of the spherical cavity in which the black hole is enclosed. 
Now note that 



w s 'Uh s )-w {E) - T "^ E) - T '^ s (58) 



One more relation to be noted is that derivable from the expression for (3, 



Using these equations we get 



d JR V(R-2M) 

dp M ~^ (R-3M) (59) 



2 1 R 

a ^^3M~R (60) 



3 mVR-2M 

Q 3 = - C/ ,_ 2 , D oTF\3 ( 61 ) 



64tt 2 (R-3M) 3 
-3 R 2 (9M-4R ) 
(8^)3 (i?-3M) a 



"4= 7^3 TWZ^mVs ( 62 ) 



Plugging in Ea. H31(l . The final Expression comes out to be 

S = S - M^)-) + + i h (A^f . 2 ) - l ^S*™ (63, 

This expression till first order was derived in |3fij . One should note here that the first order corrections are not of 
the type ln(S ). The higher order corrections are also not of the form 1/S but they have the same order as 1/Sq. 



IV. DISCUSSIONS 



In this paper we have developed a scheme to calculate the corrections to any order for a thermodynamic system. 
The higher order corrections arise when one expands the canonical entropy function around the extremum. The 
correction when applied to the BTZ black hole in the AdS Schwarzschild black hole results in terms that can be 
arranged to be a polynomial in I/Sq. For the Schwarzschild black hole in a cavity the correction higher in order than 
the logarithmic correction is of order 1/Sq. The scheme used is universally applicable as it is model independent. 

Since no quantum ideas were used anywhere in the derivation, the corrections calculated are corrections to the 
microcanonical entropy in a classical theory. Where such corrections have been found using quantum gravity theories, 
the leading order corrections are logarithmic and hence greater than the higher order corrections (1/S) found here |l6j |. 

Now we would like to bring into light some assumptions that were made while doing the calculations. 

The Bekenstein Hawking entropy is taken to be the canonical entropy and it is used while calculating the 
derivatives for getting the corrections. We expect the canonical entropy to have an extremum to which the system 
will settle to. But this is not so for the Bekenstein Hawking Entropy. We demonstrate this by a specific example. 
The Bekenstein Hawking entropy for the BTZ black hole has a 4 dependence and hence the extremum corresponds to 
j3 = oo. To the first approximation the canonical entropy and the Bekenstein Hawking Entropy agree |17| . The way 
out of this situation might be as pointed out in p"sl | inspired by Conformal Field Theory that the complete entropy 
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function may be having other terms which go to zero in the limit that the black hole has mass far greater than the 
Planck Mass; but the function then has a saddle point due to those extra terms. One such proposed function was 



S<J3) = ap m + A (64) 



In such a case the extremum occurs at 



/9b = ( -) ^ (65) 
If we proceed along the lines of 0] and expand S around /3q we get 

S(/3) = a{a n b m )T^ + I 7 ( a "+ 2 6™- 2 )^ (/3 - p Q f + 5(a n+3 b m - 3 ) ((3 - /3 ) 3 + C0" +4 & m ~ 4 ) V&V (p _ /? ) 4 + ... 



where a, 7, 5, £ are functions of m and n. 

We then proceed to use Ea. ^Hlf) to get the entropy 



(66) 



S=S - iln(5 T|)+/(m,n)A + ... ( 67 ) 
z i 

thus verifying the prediction in [is[ that the subleading order terms have dependence on m and n (order unity) 
and also that such an entropy function has a subleading order correction of the form One should note that the 
logarithmic corection also has a term involving m and n and it has been separated out. 

The analysis shall go through only under the circumstances of having such an entropy function to provide for an 
extremum and it shall breakdown when the Planck mass limit is approached. 

The logarithmic correction due to a non-uniform energy spectrum of the black hole results in a logarithmic 
correction in addition to the other corrections. Such a correction was first noted in It would be wrong to discard 
such spectral factors apriori without having a complete theory of spacetime fluctuations. We have left this factor in 
the final expression so as to allow for such conditions if they need to be considered. This factor can prove important 
in deciding whether the system is thermally stable or not as has been shown in |37| . 

Because the BTZ entropy function can be obtained from AdS Schwarzschild one by putting d — 3 one can verify 
that the correction term in Ea. H54l) is the same as in Ea. l42|) if one substitutes d — 3. It is interesting to note that 

for the BTZ blackhole, the correction that we found out is negative. The correction is also negative for the 
AdS Schwarzschild black hole for d > 3. At this moment it would be bold to conjecture the universality of this sign 
though. 

Also to be noted is the fact that the order of the neglected terms in Eq. i|31[l is dimensionally of the order i . 

This means that corrections to the order of can also be found out using Ea. (|31l) . Similarly all the higher order 



corrections can be calculated using Ea (|22[) . 

We have seen that the thermal fluctuations result in corrections to the black hole entropy. As a result the Cardy 
Vcrlinde formula receives corrections via the holographic relations. The effect of the logarithmic corrections on the CV 
formula was studied in [33,j33. One can in a similar manner find out the corrections that the CV formula acquires 
in light of the higher order corrections that have been derived in these papers. 

The holographic principle can be tested by examining the higher order corrections to black hole entropy. The 
general entropy corrections calculated by the Cardy formula were compared with the corrections to the black hole 
entropy from canonical considerations in j2]j to test the holographic principle. The correction to the black hole 
entropy derived in [27j differ for both the BTZ blackhole and the AdS blackhole. The reason for the same was pointed 
out earlier. 
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